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I. INTRODUCTION 



In the Standard Model (SM), the charged-current weak interaction hnks the top quark 
with a down-type quark with an amphtude proportional to the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix Vtq {q = d, s, b). At hadron colliders, this interaction leads to a single top 
quark final state through three possible processes: the s-channel decay of a virtual W 
{qq' W* — > tb), the t-channel exchange of a virtual W {bq tq' and bq' tq, also 
referred to as VT-gluon fusion), and the associated production of a top quark and a W 
boson [bg tW~). Since the cross section for each of these processes is proportional to 
the CKM matrix element \Vtb\^, a measurement of each single top quark production cross 
section determines \Vtb\- A study of spin correlations in single top quark production can be 
used to test the V — A structure of the top quark charged-current weak interaction. Such 
structure implies that the top quark should be polarized. 

Among the three processes, the largest cross section at both the Tevatron proton- 
antiproton collider at Fermilab and the proton-proton Large Hadron Collider (LHC) at 
CERN is due to the t-channel. The next largest cross section is from the s-channel at 
the Tevatron, and from associated production at the LHC. The s-channel production cross 
section is small at the LHC because it involves a quark- ant iquark collision. Similarly, as- 
sociated production cross section is relatively large at the LHC because the gluon parton 
distribution functions (PDF) grow more rapidly with decreasing x than the light quark PDF. 
The s-channel and t-channel processes should be observed for the first time in Run II at 
the Tevatron; whereas the observation of associated production will likely have to await the 
LHC. 

Single top quark production is also a very important background to many searches for new 
physics. For example, the s-channel process is a significant background to Hi^s searches 
at the Tevatron in the production process qq' —> WH with decay if — > 66 and 
other new physics searches j^. At the LHC, the Wt associated jproduction is an important 
background to Higgs searches „r the decay channel H WwU and the primary charged 
Higgs boson production channel bg — > tH^ with ^ tz/ p|. As the largest single top 
process, the t-channel production is an important background to many new physics searches. 

In Tevatron Run I, searches for single top quark production were performed by both the 
D0 0] and CDF collaborations j^. At the 95% confidence level, the D0 limit on the s- 
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channel production is 17 picobarn (pb) and the CDF hmit is 18 pb. At the same confidence 
level, the D0 limit on the t-channel production cross section is 22 pb and the CDF limit 
is 13 pb. Searches for single top quark production in Run II have begun, and the limits 
from CDF are 13.6 pb in the s-channel and 10.1 pb in the t-channel The s-channel and 
t-channel single top quark processes can be probed separately at the Tevatron by taking 
advantage of 6-quark tagging using displaced vertices and differences in the kinematic distri- 
butions of the fe-tagged and non-6-tagged jets. Usually, only one 6-tagged jet can be expected 
in the t-channel case while two 6-tagged jets can be expected in the s-channel case. This is 
because the b quark produced with the top quark tends to be collinear with the initial state 
gluon in the t-channel, giving it a large pseudo-rapidity (77) and low transverse momentum 
[Pt) and thus making it challenging to 6-tag this jet experimentally. It is important to 
separate the two processes since they have different systematic uncertainties and different 
sensitivities to new physics |lo| . 

The extraction of a signal is more challenging for single top quark production than for 
top pair production since there are fewer objects in the final state and the overall event 
properties are less distinct from the large VT+jets background. Therefore, an accurate 
calculation including higher order quantum chromo dynamics (QCD) corrections is needed. 
The next-to-leading order (NLO) Ofa J corrections to single top quark production have 



already been carried out in Refs. (ll|. 



12| . which shows an uncertainty on the total cross 



section of about 5% by varying the factorization and renormalization scales. For the s- 
channel, multiple soft gluon resummation effects have been calculated In order to 

confront theory with experimental data, where kinematical cuts are necessary in order to 
detect a signal, it is crucial to accurately model event topologies of single top quark events. 

cross sections for on-shell single top 
quark production. The complete NLO calculations including both the single top qu ark 
production and decay have been done independently by two groups recently [10, ll2l ll^l • In 
both calculations, the narrow width approximation was adopted in order to link top quark 
production with its consequent decay 3, Q| • In Ref . Q] , various kinematics distributions 
are examined both with and without top quark decay at NLO. In our previous study |l8|. we 
presented a detailed phenomenological analysis of s-channel single top quark production at 
the Tevatron, focusing on signal cross sections and kinematical distributions. To complete 
our study of single top processes at the Tevatron, we will focus on the single top quark 
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production process with the largest production cross section in this paper, namely the t- 
channel. A more realistic study which disentangles the signal of single top quark events from 
copious backgrounds and further separates the s-channel from the t-channel will be given in 
the future. 

The paper is organized as follows. In SeclHl we first present the inclusive cross section for 
the t-channel single top process using both the fixed mass narrow width approximation and 
the "modified" narrow width approximation. We then investigate its dependence on the top 
quark mass as well as the renormalization and factorization scales. In Sec. IIIIl we examine 
the acceptance of single top quark events for various sets of kinematic cuts. In Sec. IIV| 
we explore the kinematical information in the final state objects. Section contains our 
conclusions. 



II. CROSS SECTION (INCLUSIVE RATE) 

In order to accurately measure Vtb using the single top process, the uncertainties in the 
theoretical prediction need to be minimized. In this section, we show the inclusive production 
rates and discuss their dependence on the factorization and renormalization scales and the 
top quark mass. 

We present numerical results for the production of single top quark events considering 
the leptonic decay of the VT-boson from the top quark decay at the upgraded Tevatron 
(a 1.96 TeV pp collider). Unless otherwise specified, we use the NLO parton distribution 
function (PDF) set CTEQ6M Q], defined in the MS scheme, and the NLO (2-loop) running 
couphng as with provided by the PDFs. For the CTEQ6M PDFs, A^^ = 0.326 GeV 
for four active quark flavors. The values of the relevant electroweak parameters are: a = 
1/137.0359895, G„ = 1.16637 x 10"^ GeV-^ Mw = 80.40 GeV, Mz = 91.1867 GeV, and 



sin^6'iy = 0.231 j20|. The square of the weak gauge coupling is g'^ = 4a/2Mm/Gu- As i 



m 



the s-channel study here we also focus our attention only on the positively charged 
electron (i.e., positron), though our analysis procedure also applies to the muon final state. 
Including 0{as) corrections to W ^ qq', the decay branching ratio of the W boson into 
leptons is Br{W — > I'^v) = 0.10 8 | 20j] . Unless otherwise specified, we will choose the top 
quark mass to be = 178 GeV |2ll . l22l | and the renormalization scale (yU/j) as well as the 
factorization scale (fJ^p) to be equal to the top quark mass. The top quark mass and scale 
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dependencies of the t-channel production cross section are investigated in this section. 



A. Theoretical CutofF Dependence 

The NLO QCD differential cross sections are calculated using the one-cutoff phase space 
slicing (PSS) method This procedure introduces a theoretical cutoff parameter 

(■5mm) in order to isolate soft and collinear singularities associated with real gluon emission 
subprocesses by partitioning the phase space into soft, collinear and hard regions such that 

In the soft and collinear regions the cross section is proportional to the Born-level cross sec- 
tion. Using dimensional regularization, we can evaluate the real gluon emission diagrams in 
n-dimensions under the soft gluon approximation in the soft region, or the collinear approx- 
imation in the collinear region, and can integrate out the corresponding phase space volume 
analytically. The resulting divergences are cancelled by virtual corrections or absorbed into 
the perturbative parton distribution functions in the factorization procedure. Since the cut- 
off is introduced in the calculation only for this technical reason and is unrelated to any 
physical quantity, the inclusive rate must not depend on it. In other words, the sum of all 
contributions, virtual, resolved, and unresolved corrections must be independent of Smin- 
This is the case as long as Smin is small enough so that the soft and collinear approximations 
are valid. However, numerical cancellation in the Monte Carlo (MC) integration becomes 
unstable if Smin is too small. Furthermore, the jet-finding algorithm and other infrared-safe 
experimental observables should also be defined in a way such that they are consistent with 
the choice of Smin- In practice, one wants to choose the largest Smin possible within these 
constraints in order to minimize the processing time of the MC integration program. For 
our study, we found a value of Smin = 1 GeV^ to be appropriate for studying the t-channel 
single top process. For comparison, in Ref. Q|, we found that a value of Smm = 5 GeV^ is 
adequate for studying the s-channel single top process. A detailed discussion of the phase 
space slicing method can be found in Ref. [Hf. 

Figure ^ illustrates that the total t-channel production cross section is indeed insensitive 
to Smin for a large range of Smin values. The figure shows the sum of the virtual and 
unresolved real corrections (s + v) as well as the resolved contribution (real) to the t-channel 
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FIG. 1: The theoretical cutoff Smin dependence of the inclusive t-channel single top quark cross 
section at the Tevatron with ^ur = fip = nit for uit = 178 GeV. The decay branching ratio 
t — > bW^i^ e^u) has been included. 



single top process as a function of Smin. Although contributions from the individual pieces 
vary, their sum (total) remains essentially constant for a large range of Smin- With the choices 
of = /ii? = at the Tevatron, we obtain an inclusive cross section for the t-channel 
single top (t only) processes (with M^-boson decay branching ratio) of 104.8 femtobarn (fb), 
which agrees with Ref. Q]. 



B. Inclusive Cross Section 



To facilitate the calculation, we divided the higher-order QCD corrections into three 
separate gauge invariant sets: corrections from the light quark line of top quark production 
(LIGHT), corrections from the heavy quark line of top quark production (HEAVY), and 
corrections from the top quark decay (TDEC). The explicit diagrams and definitions for 
these three sets can be found in Ref. [l7|. The inclusive cross section as well as the individual 
0{as) contributions are listed in Table HI The effects of the finite widths of the top quark 
and H^-boson have been included. As one expects, the dominant 0{as) corrections come 
from the heavy quark line, because the bulk part of the radiative corrections originating 
from the light quark line have been absorbed into the definition of the light quark PDFs. 
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Cross Section 


Fraction of 




(fb) 


NLO (%) 


Born -level 


99.2 


94.6 


o(q,) heavy 


5.56 


5.31 


0{as) LIGHT 


1.03 


0.98 


0(a,) TDEC 


-0.81 


-0.77 


0{as) sum 


5.54 


5.28 


NLO 


104.8 





TABLE I: Inclusive t-channel single top production cross section for different subprocesses, includ- 
ing the top quark decay branching ratio t bW~^{-^ e'^v). 

For the same reason pointed out in our previous paper the TDEC contribution is very 
small compared to the other two. 



C. Top Quark Mass and Renormalization/Factorization Scale Dependence 

In order to test the SM and measure the CKM matrix element Vtb-, one needs an ac- 
curate prediction of single top quark production and decay so as to reduce the theoretical 
uncertainties. Besides the top quark mass, the choices of renormalization and factorization 
scales also contribute to the uncertainty of the theoretical prediction. The renormalization 
scale yUij is introduced when redefining the bare parameters in terms of the renormalized 
parameters, while the factorization scale /i^ is introduced when absorbing the collinear di- 
vergences into parton distribution functions. Therefore, both /i^ and are unphysical and 
the final predictions should not depend on them. However, since we work at a fixed order 
in perturbation theory, we indeed see a dependence of the predicted cross section on /i/j and 
/ii?. The change due to varying the scale is formally of higher order. Since the single top 
quark rate is small at the Tevatron, it is important to reduce the scale uncertainty in order 
to compare the theory prediction with experimental data. Here, we examine the top quark 
mass and the scale dependence of the t-channel production cross section. 

As shown in Fig.|2l the cross section changes by about ±9% when the top quark mass irtt 
is varied by its current uncertainty of about =i=5 GeV around 178 GeV. It can also be seen 
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FIG. 2: Top quark mass dependence as well as renormalization and factorization scale dependence 
of the inclusive t-channel single top quark cross section at the Tevatron. The decay branching ratio 
of t — > bW~^{-^ e'^u) has been included. 

that measuring the top quark mass to an uncertainty of 1 — 2 GeV will reduce the theoretical 
uncertainty on the single top cross section correspondingly. The reduction in uncertainty on 
the t-channel single top production rate will improve the measurement of the CKM matrix 

In order to examine the scale dependence of the t-channel single top production rate, we 
show in Fig.|21the results of two typical scales: one is the top quark mass {iip = I^^r"^ = 'm^t)-, 
shown as the solid line, the other is the W boson mass {^p = ^^^B°'^ = My/), shown as the 
dashed line. For the decay of top quark, we take /i^^*^ = m^, which gives similar results as 
the choice of ^'^^'^ = My/. The band constrained by these two scales represents a range 
of uncertainty due to the NLO predictions. The usual practice for estimating the yet-to-be 
calculated higher order QCD correction to a perturbative cross section is to vary around the 
typical scale by a factor of 2, though in principle the "best" scale to be used for estimation 
cannot be determined without completing the higher order calculation. In Fig. El we show 
the total cross section of the t-channel single top production for a range of the scale 
which, for simplicity, is set to be equal to ^r. We have multiplied the LO cross sections by 
a constant factor (the "K-factor") of 1.04 (shown as the dashed line) in order to compare 
to the NLO ones (shown as the solid-line). It is clear that the NLO calculation reduces 
the scale dependence. For example, when the scale is changed from Hf/'^ to 2jjip around 
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FIG. 3: Inclusive t-channel single top quark production cross section at the Tevatron for mt = 178 
GeV, versus the ratio of the factorization scale fip to its typical value fiQ, where ^uq = mt (solid 
line) and /xq = M\y (dashed line), respectively. The decay branching ratio of t — > bW~^{^ e^u) 
has been included. 



= f^t, the LO rate varies by —0.9% to —1.0% while the NLO rate varies only by —0.6% 
to +1.3%. Similar results also hold for varying the scale around hf = Mw- 



III. SINGLE TOP ACCEPTANCE STUDIES 



In this section we explore the final state objects of t-channel single top quark events. 
Although the IV-boson from the top quark decay can decay in both leptonic and hadronic 
modes, we focus only on the leptonic decay mode in this study because the all-jet produc- 
tion mode of single top quark events is difficult to observe experimentally. Therefore, the 
signature of t-channel single top quark events which is accessible experimentally consists of 
one charged lepton, missing transverse energy, together with two or three jets. Since we 
are studying the effects of NLO QCD radiative corrections on the production rate and the 
kinematic distributions of single top quark events at the parton level in this paper, we do 
not include any detector effects, such as jet energy resolution or 6-tagging efficiency. Only 
an approximation of kinematic acceptances of a generic detector are considered. Since the 
Tevatron is pp collider, and pp is a CP-even state, the production cross sections for bi{j) at 
the Tevatron are the same as those for tb{j) (when ignoring the small CP violation effect 
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FIG. 4: Pictorial illustration of the notation used in this paper. 



induced by the CKM mixing matrix in quarks). For simplicity we therefore only consider 
distributions of t-channel single top quark events in this work in which a t quark (not in- 
cluding i) decays into W^{-^ l^v) and a h quark. In this section, we first present the event 
topology of the t-channel single top process, and then introduce a jet finding algorithm and 
the various kinematical cuts in order to study the acceptance. 

A. Event Topology 

At tree level, the collider signature of the t-channel single top process includes two jets 
(one 6-tagged jet from the h quark from the top quark decay, and one non-6-tagged jet from 
the light quark), one charged lepton, and missing transverse energy {j^t) in the final state. 
This signature becomes more complicated beyond tree level, as Fig. |3] indicates. The light 
quark jet is also called "spectator jet", and the label "untagged jet" refers to all jets which 
do not contain a 6 or 6 quark. 

At NLO, besides the charged lepton and ^t, there may be two jets (one 6-tagged jet and 
one untagged jet) as for the Born-level, or there may be three jets. The flavor composition 
of the three-jet flnal state depends on the origin of the third jet. When it is a gluon or 
anti-quark (cases (a-c) in Fig. there is one 6-tagged jet and two untagged jets. When 
it is a 6 quark (case (d) in Fig. El also called H^-gluon fusion), there are two 6-tagged jets 
and one untagged jet. Therefore, prescriptions are needed to identify the h jet from the top 
quark decay and the light quark jet produced with the top quark. 

In our study, we differentiate the following three cases: 
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FIG. 5: Representative diagrams of the real emission corrections for the t-channel single top process: 
(a) and (b) represent the real radiative corrections to the LIGHT quark line, while (c) and (d) 
represent the real radiative corrections to the HEAVY quark line. The NLO QCD corrections are 
indicated by the large shaded ellipse. Detailed Feynman diagrams can be found in Ref. [l^ . 

1. Born- level-type exclusive two-jet events (containing the b quark and the light quark): 
In Born-level type events, i.e. without an additional hard jet, the final state is given 
by the b jet from the top quark decay and the light quark jet produced with the top 
quark. Experimentally, this configuration is easily reconstructed because the 6-tagged 
jet is identified as the b quark, and the untagged jet is assigned to the spectator jet. We 
note that in our phenomenology study to be presented in Sec. II VI we shall concentrate 
on inclusive two-jet events and exclusive three-jet events. 

2. Exclusive three-jet events with one b jet (containing the b quark, the light quark, and 
a gluon or anti-quark): 

When the radiated gluon is reconstructed as a separate jet in the LIGHT, TDEC, or 
HEAVY corrections, the final state contains one b jet and two other jets. As shown 
later, the transverse momentum and energy differences can be used to separate the 
spectator jet from the gluon jet (third jet). The transverse momentum of the spectator 
jet, which comes from the initial quark (g) after emitting the effective VT-boson, peaks 
around ~ Mvi//2, thus its energy is large. In comparison, the transverse momentum of 
the gluon jet is small. Due to the coUinear enhancement, the resolved gluon jet prefers 
to move along the beam line direction at smaller transverse momentum {px) for both 
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the LIGHT and HEAVY corrections, and along the b quark moving direction for the 
TDEC correction. 



3. Exclusive three-jet events with one b jet and one b jet (containing the b quark, the 
light quark, and a b quark): 

This final state is produced in the ly-gluon fusion process, case (d) of Fig. El In 
this case, the spectator jet can be uniquely identified, but the b jet from the top 
decay cannot because heavy fiavor tagging algorithms do not distinguish b jets from 
b jets experimentally. Although the likelihood of tagging two 6-jets is smaller than 
that of tagging a single fc-jet, requiring two 6-tagged jets will suppress the QCD and 
VT+jets backgrounds significantly. We again can use transverse momentum differences 
to separate b jets from b jets. The b jet from the top quark decay has a transverse 
momentum peaking around mj/3, while the b jet from gluon splitting tends to move 
along the beam line (the gluon-moving direction) due to the collinear enhancement and 
is much softer. These kinematic differences also enable us to separate the t-channel 
process from the s-channel process, recalling that in the s-channel process, both the 
b jet and the b jet are preferentially produced at central rapidity and large pt- 

The unique signature of the t-channel single top process is the spectator jet in the forward 
direction, which can help to suppress the copious backgrounds, such as Wbb and ti produc- 
tion. Studying the kinematics of this spectator jet is important in order to have a better 
prediction of the acceptance of t-channel single top quark events and of the distribution of 
several important kinematic variables. In this work, we study the impact of the NLOQCD 



corrections on the kinematic properties of the spectator jet. As pointed out in Ref. [^j, m 
the effective-iy approximation, a high-energy t-channel single top quark event is dominated 
by a longitudinal W boson and the b quark fusion diagram. It is the same effective lon- 
gitudinal W boson that dominates the production of a heavy Higgs boson at high energy 
colliders via the ly-boson fusion process. For a heavy SM Higgs boson, the longitudinal 
W boson fusion process dominates the Higgs boson production rate. Therefore, it is also 
important to study the kinematics of the spectator jet in t-channel single top quark events 
in order to have a better prediction for the kinematics of Higgs boson events via the WW 
fusion process, i.e. qq(WW) Hq'q' at hadron colliders. 



12 



B. Acceptance 



In order to meaningfully discuss the effects of gluon radiation in single top quark events, 
we must define a jet as an infrared-safe observable. In this study, we adopt the cone-jet 
algorithm j^l as explained in our previous work More specifically, we adopt the E- 
scheme cone-jet approach (4-momenta of particles in a cone are simply added to form a 



jet) with radius R = yArp + A^^ in order to define 6, q and possibly extra g, q, or b jets, 
where Arj and A0 are the separation of particles in the pseudo-rapidity rj and the azimuthal 
angle 0, respectively. For reference, we shall consider both R = 0.5 and R = 1.0. The same 
i?-separation will also be applied to the separation between the lepton and each jet. 
The kinematic cuts imposed on the final state objects are: 

> 15GeV , \7]e\ < r/f"^, 

> 15 GeV , 

> 15 GeV , \rij\ < r/f'^", 

ARij > Rcut , ^Rjj > Rcut, (2) 

where the jet cuts are applied to both the b- and light quark jets as well as any gluon or 
antiquark jet in the final state. Two lepton pseudo-rapidity cuts are considered here: a loose 
version with -q^"-^ = 2.5 and a tight version with jj^"-^ = 1.0. Similarly, loose and tight cuts 
are also considered for the jet pseudo-rapidity, r^J*"^ = 3.0 and r^™"^ = 2.0, respectively. The 
minimum transverse energy cuts on the lepton as well as the jets is 15 GeV. Each event is 
furthermore required to have at least one lepton and two jets passing all selection criteria. 
The cut on the separation in R between the lepton and the jets as well as between different 
jets is given by Rcut which is is chosen to be 0.5 (or 1.0). 

In Table [HI we show the single top production cross sections in fb for the loose and tight 
set of cuts listed in Eq. (0) for the different subprocesses. A larger value for R^ut reduces 
the acceptance significantly mainly because more events fail the lepton-jet separation cut. 
While this is only a 13% reduction at the Born-level, the difference grows to 25% at NLO. 
Hence, a smaller cone size is preferred in order to keep the acceptance at a high level. For 
events with at least three jets, imposing a harder cut on the transverse momentum of each 
jet {Etj > 30 GeV) decreases the contribution from the LIGHT and HEAVY corrections by 
a factor of 3 to 4, but the contribution from the TDEC correction by more than a factor 
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a[fb] 


LO 


NLO 


Heavy 


Light 


Decay 


Tevatron (a) 


tq + tqj 


65.6 


64.0 


4.9 


-3.4 


-0.59 




tq + tqj, Exi > 30GeV 


46.3 


44.0 


4.2 


-3.3 


-1.7 








25.3 


15.5 


5.1 


4.8 




tqj, En > 30GeV 




7.4 


5.4 


1.3 


0.75 


Tevatron (b) 




56.8 


48.1 


-0.45 


-4.4 


-2.2 




tq + tqj, Exi > 30GeV 


40.1 


33.3 


0.0 


-3.3 


-2.1 








14.8 


10.8 


2.2 


1.7 




tqj, En > 30 GeV 




4.5 


3.7 


0.58 


0.20 


Tevatron (c) 


tq + tqj 


31.1 


34.0 


5.8 


-2.8 


0.82 




tq + tqj, Etj > 30 GeV 


24.4 


24.2 


3.8 


-2.7 


-0.7 








11.4 


6.7 


2.4 


2.3 




tqj, Etj > 30 GeV 




3.9 


2.8 


0.74 


0.41 



TABLE II: The f-channel single top production cross section at the Tevatron (in fb ) for different 
subprocesses under various cut scenarios: (a) is the loose version with rj'P'^^ = 2.5, rjj^'^^ = 3.0, and 
Rcut = 0.5, (b) is the loose version with a larger jet clustering cone size (and jet-lepton separation 
cut) of Rcut = 1-0, and (c) is the tight version of cuts with ry™"^ = 1.0, 7]"^°-^ = 2.0, and Rcut = 0.5. 
The first two columns show the Born-level and full NLO cross sections, the last three columns show 
the contributions from the different 0{as) processes. The decay branching ratio t 6VF(— > ev) is 
included. 

6. This large reduction occurs because the top quark mass sets the scale for the top quark 
decay contribution rather than the invariant mass of tq' system, resulting in a softer jet Et 
spectrum for the TDEC correction (cf. Fig. |211 to be discussed in Sec. IIVD|) . 

Figure IHl shows how the observed cross section chang GS clS db function of the jet Et cut 
when applying the loose set of cuts, including a requirement of there being at least two jets 
in the event. The figure also shows the dependence of the fraction of two-jet events and 
three-jet events on the jet Et cut. At the Born-level, there are only two-jet events, whereas 
0{as) corrections can produce an additional soft jet. The fraction of events with these 
additional jets is low only for very high jet Et thresholds. For typical jet Et thresholds 
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jet Pj threshold [GeV] jet p,. threshold [GeV] 

(a) (b) 

FIG. 6: Cross section and fraction of three-jet events at NLO for varying jet pj- cuts, requiring 
only that rijets ^ 2, and not making any cuts on the electron or neutrino. Shown is the total 
cross section for events with two or three jets as a function of the jet Ex cut for three different jet 
pseudo-rapidity cuts (a) and the fraction of three-jet events as a function of the jet px for different 
jet pseudo-rapidity cuts (b). The lowest threshold considered is 15 GeV. 

considered by experiments of 15 GeV to 25 GeV these jets add a significant contribution. 
As expected, the effect is not quite as large when only jets within a very small t] range are 
considered because the extra jet typically has higher 77. In order to study 0{as) effects it is 
thus important to set the jet rj cut as high as possible and the jet Ej- cut as low as possible. 

As mentioned before, the event rate for single top quark events is small and it is important 
for experiments to maximize their acceptance. We will thus use the loose set of cut values 
for the following discussion: r;™'^^ = 2.5, 77™'^^ = 3.0, and Rcut = 0.5, E^J"- = 15 GeV, cf. 
Eq. (0). 



IV. SINGLE TOP QUARK EVENT DISTRIBUTIONS 

In this section we examine the kinematic properties of single top quark events. As dis- 
cussed in the previous section, the signature of the t-channel process includes at least one 
6-tagged jet, one untagged jet, one charged lepton, and missing energy. At the Born-level, 
it is straightforward to identify the single top final state using the 6-tagged jet and the 
reconstructed W boson to reconstruct the top quark and identifying the spectator jet as 
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the light quark. At NLO, however, an additional jet is radiated, which will complicate the 
reconstruction of the top quark final state. First, this is because the additional jet can be 
either a 6-tagged jet or an untagged jet. When it is the 6-tagged jet, we need to select 
the correct b quark from the two possible jets in the finial state to reconstruct the top 
quark. When it is the untagged jet, we need to select the correct spectator jet. Second, 
the additional untagged jet can come from either the production or the decay of the top 
quark. Production-stage emission occurs before the top quark goes on shell and decay-stage 
emission occurs only after the top quark goes on shell. In production emission events, the 
W boson and b quark momenta will combine to give the top quark momentum, while in the 
decay emission event the gluon momentum must also be included in order to reconstruct 
the top quark momentum properly. To find the best prescription for identifying the correct 
b jet and spectator jet, we first examine various kinematic distributions of the final state 
particles. We then investigate two top quark reconstruction prescriptions: using the leading 
6-tagged jet and the best jet algorithm. Choosing the 6-tagged jet prescription, we can also 
improve the reconstruction efficiency for the W boson. We then study the effects of NLO 
corrections on distributions concerning the reconstructed top quark, in particular spin cor- 
relations between the final state particles. Finally, we explore the impact of the radiated jet 
in exclusive three-jets events. We use only the loose set of cuts to maximize the efficiency 
when examining the distributions and efficiencies in detail. 

A. Final State Object Distributions 

1. Charged Lepton and Missing Transverse Energy 

In this section we examine various kinematic distributions of final state objects after 
event reconstruction and after applying the loose set of cuts, cf. Table H] and Eq. (jSj). 
We concentrate on inclusive two-jet events in this section because they give more reliable 
infrared-safe predictions. 

Figures [7| (a) and (c) show the transverse momentum of the electron and the missing 
transverse energy ^t, respectively. As expected, because they are leptons and not quarks, 
the change in shape when going from Born-level to NLO is not very large for the electron or 
the missing transverse energy. The pseudo-rapidity distribution of the electron is given in 
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FIG. 7: The transverse momentum pT (a) and pseudo-rapidity t] (b) distributions of the electron 
and the missing transverse energy (c) after selection cuts, comparing Born-level to 0{as) 
corrections. 

Fig-Etb). This distribution widens at NLO, but again the effect is small because the overall 
0{as) contribution to the event rate is small, cf. Table H] We note that the peak position 
of the distribution is at a higher value than that of the electron because the neutrino 
from the IV-boson decay moves preferentially along the direction of the top quark. This 
is due to the left-handed nature of the charged current interaction and can easily be seen 
when examining the spin correlations between the charged lepton and the top quark in the 
top quark rest frame. We will comment more on this subject in Sec. II V (^1 



2. Spectator Jet 

The differences between the s-channel and t-channel single top processes at the Born-level 
are that the former has two 6-taggable jets in the final state (the b quark and the b quark) 
while the latter only has one 6-taggable jet and also has one light quark jet. This untagged jet 
(spectator jet) is a unique feature of the t-channel process which can be used to disentangle 
t-channel single top quark events from the copious backgrounds. Therefore, its kinematic 
distributions need to be well studied, especially the impact of 0{as) corrections on the 
features which make this jet a unique signature, such as the spectator jet pseudo-rapidity, 
cf. Fig. El 

The pseudo-rapidity distribution of the spectator jet is asymmetric because the Tevatron 
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FIG. 8: Transverse momentum pT (a), energy E (b) and the pseudo-rapidity r] (c) of the spectator 
jet after selection cuts, comparing Born-level to 0{as) corrections. 



is a pp collider j25|. In order to produce a heavy top quark decaying to a positively charged 
lepton, the valence quark from the proton is most important, implying that the light quark 
will tend to move in the proton direction. We define the positive z-direction to be the 
proton direction in the laboratory frame, thus the pseudo-rapidity of the spectator jet will 
tend be positive. Similarly, the spectator jet in an anti-top quark event produced from the 
t-channel process will preferably be at a negative pseudo-rapidity due to the large anti-up 
quark parton distribution inside the antiproton. The 0{as) corrections shift the spectator 
jet to even more forward pseudo-rapidities due to additional gluon radiation. However, since 
the 0{as) corrections are small compared to the Born-level contribution, the spectator jet 
pseudo-rapidity distribution only shifts slightly. As Fig. El shows, the LIGHT and HEAVY 
contributions have almost opposite behavior. The former shifts the spectator jet to even 
higher pseudo-rapidities, while the later shifts it more to the central rapidity region. This 
behavior is due to two different effects, as illustrated in Fig. (b), in which "PA" denotes 
that the light quarks come from the proton while the bottom quarks from the anti-proton 
and vice versa for "AP". After separating the contributions by whether the light quark is 
from the proton or the antiproton, it can be seen that the HEAVY corrections shift the 
proton contribution down and the antiproton contribution up due to the slight change in 
acceptance caused by the additional jet. The LIGHT corrections show the opposite tendency. 
For the TDEC contribution, all corrections have similar shapes and the sum of them leaves 
the spectator jet pseudo-rapidity unchanged, as expected. After summing the negative soft- 



18 



-4 



LIGHT 



-1 



TDEC 



-2 







2 



-2 







2 



Jet. 



spectator 



Jet. 



spectator 



(a) 



(b) 



FIG. 9: Each individual contribution of the 0{as) corrections to the spectator jet pseudo-rapidity, 
summed (a), separately for the case when the incoming up-type quark is from the proton and anti- 
proton (b). Here, "PA" and "AP" denotes the initial state light quark originating from proton and 
anti-proton, respectively. 

plus-virtual corrections with the real emission corrections, we obtain the result shown in 
Fig. IHl which shows that the 0{as) correction shifts the spectator jet to be in even more 
forward direction. 

Besides its forward rapidity, the spectator jet also has large transverse momentum. Since 
it comes from the initial state quark after emitting the effective W boson, the transverse 
momentum peaks around ~ Mw/2, cf. Fig. |H| By comparison, the third jet is most often 
much softer, we can thus use pt of the jet to identify the the spectator jet when considering 
exclusive three-jet events. 



Compared to the lepton and ^t, the effects of the 0{as) corrections on the reconstructed 
b jet are more pronounced. Figure ITUl shows a comparison of the b jet Pt distribution between 
the Born-level and 0{as) corrections. The pt distribution of the b jet is predominantly de- 
termined by the top quark mass and therefore peaks at ~ mt/S. The NLO QCD corrections 
broaden the transverse momentum distribution and shift the peak position to lower values. 
The location of the mean of the distributions depends on the cuts that are applied because 
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FIG. 10: Transverse momentum px of the h jet after selection cuts, comparing Born-level to 0{aa) 
corrections (a) and the individual 0{as) contributions (b). 




FIG. 11: Pseudo-rapidity r] of the 6-jet after selection cuts, comparing Born level to 0{as) correc- 
tions. 

the different 0{as) corrections have different effects: in the case of the loose cuts it increases 
from 64 GeV at the Born-level to 62 GeV at full NLO. The LIGHT corrections shift the 
mean of the h jet pt distribution up, the HEAVY corrections leave it mostly unchanged, 
and the TDEC corrections tend to shift it down. The h jet px distribution receives a large 
contribution from the TDEC corrections, as expected. When a gluon is radiated from the 
top quark decay, it tends to move along the h jet direction due to coUinear enhancements 
and therefore shifts the h jet pt distribution to the small pt region, as shown in Fig. ITUl 

The b jet pseudo-rapidity distribution is less affected by the 0{as) corrections, as can be 
seen in Fig. ^2 The top quark is so heavy that it is mostly produced in the central rapidity 
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FIG. 12: Momentum of the b and b jets (a) and fraction of events with one or two 6-tagged jets as 
a function of the jet pT threshold, for both inclusive two-jet and exclusive three-jet events (b). 

region and thus the b jet from its decay also peaks around a pseudo-rapidity of zero. The 
shape of the b jet pseudo-rapidity distribution remains almost unchanged compared to the 
Born-level because it comes from the top quark decay. 

Figure El compares the momentum of the b jet and the b jet at NLO and examines the 
fraction of events containing one or two 6-tagged jets in the final state. Events with both 
b and b jets originate from the IV-gluon fusion subprocess, qg q'bt{-^ bW{-^ ^^^))- While 
the fraction of events with two 6-tagged jets is high for low jet px-, it drops quickly. For a 
jet pt above 50 GeV, the fraction of events with both b and b jets is less than 5%. Figure IT^ 
also shows that for exclusive three-jet events, the fraction of events with two 6-tagged jets 
is much higher. About half of the events contain two 6-tagged jets for the lowest threshold 
of 15 GeV, and the fraction only goes down to 5% when the jet pt reaches about 140 GeV. 



4-. Ht Distribution 

The impact that different 0{as) corrections have on the pr of the jets is also reflected in 
event-wide energy variables such as the total transverse energy in the event {Ht), defined 
as 

Ht = p'^'"°''+:^t + T.Pt'- (3) 

jets 

The distribution of for t-channel single top quark events is shown in Fig. The Born- 
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FIG. 13: Total event transverse energy Ht after selection cuts, comparing Born level to 0{as) 
corrections (a) and the individual 0{as) contributions (b). 

level Ht distribution peaks around 200 GeV. Both the HEAVY and LIGHT contributions 
decrease the height of the peak and shift it to higher values, and all three 0{as) contributions 
broaden the distribution. 

B. Event Reconstruction 

When analyzing single top quark events we would like to take advantage not only of 
simple single-object kinematics but also of correlations between objects. In order to take 
full advantage of the correlations, we need to reconstruct the event completely, not just the 
final state jets but also intermediate particles, in particular the VT-boson and the top quark. 

The W boson can be reconstructed from the final state electron and the observed missing 
transverse energy, ]^t- The lack of information about the beam-direction component of the 
neutrino momentum (p^) that would prevent this reconstruction is typically overcome by 
requiring that the invariant mass of the electron-neutrino system be equal to the mass of the 
W boson. This additional constraint results in two possible solutions for p^. One usually 
follows the prescription given in Ref. j^l of choosing the solution which has the smaller 
This picks the correct in about 70% of the events. The price paid for this method is that 
the W boson will blur the spin correlation in the reconstructed single top quark event. We 
will see later that we can improve on this by using a top mass constraint. 
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(a) (b) 

FIG. 14: Transverse momentum of the b quark jet, for all events (solid histogram), only for those 
events for which the leading 6-tagged jet corresponds to the b quark jet (dashed histogram), only for 
those events in which the best jet corresponds to the b quark jet (dotted histogram), and only for 
those events in which the leading jet corresponds to the b quark jet (dash-dotted histogram). The 
transverse momentum is shown for the inclusive two-jet sample (a) and for the exclusive three-jet 
sample (b). 

In order to reconstruct the top quark, the reconstructed W boson then needs to be 
combined with the b jet from the top quark decay. The challenge to overcome here is the 
proper identification of the b jet. In the s-channel single top process, we used the "best-jet" 
algorithm to find the correct b jet among the two possible 6-tagged jets in the final state, 
making use of the known top quark mass (m^ = 178 GeV) [18|. The effectiveness of the 
best-jet algorithm is mainly limited by the efficiency of the W boson identification method; 
if the W boson is not reconstructed properly, then identifying the b jet from the top quark 
decay becomes a random choice. 

Unlike the s-channel processes, inclusive two-jet events of the t-channel single top pro- 
cesses typically only contain one 6-tagged jet in the final state, corresponding to the b quark 
from the top quark decay. Figure El shows a comparison of three different algorithms to 
identify the b jet from the top quark decay: a) using the leading jet (highest pt) in the event, 
b) using the leading 6-tagged jet in the event, and c) using the best jet. The leading jet in 
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this case could be either the b jet or the hght quark jet, thus the efficiency for identifying 
the fe-jet from top quark decay with method a) is rather low. The leading 6-tagged jet is 
either the b jet from the top quark decay or the b jet originating from the W^-gluon fusion 
subprocess. Figure El shows that the leading 6-tagged jet corresponds to the b quark from 
the top decay most of the time; it identifies the b quark correctly in 95% of the events in 
the inclusive sample after the loose selection cuts. We expect the efficiency to drop in the 
three-jet sample because of the additional b quark in the final state. However, as shown in 
Fig. E[ the b jet from the top quark decay is much harder than the b jet radiated from the 
initial state gluon. As a result, the leading 6-tagged jet corresponds to the b quark from the 
top quark decay in 90% of the exclusive three-jet events. Figure IT^ shows that this fraction 
is slightly smaller for low b jet px values, and reaches 100% above about 100 GeV. By com- 
parison, the best jet corresponds to the b quark in only 80% of the events in the inclusive 
sample and about 72% of the events in the exclusive three-jet sample. We therefore choose 
the method of identifying the leading 6-tagged jet as the b quark from the top decay. Having 
identified the b quark in this manner with no additional assumptions allows us to use it to 
reconstruct the top quark and the W boson more accurately, in particular to determine the 
longitudinal momentum of the neutrino (p^). We reconstruct the W boson as above using 
a W mass constraint. We then identify the proper neutrino p'^ by choosing the solution for 
which the invariant mass of the reconstructed W and the leading 6-tagged jet is closest to 
the true top quark mass, i.e. 178 GeV. 

This method identifies the correct neutrino p'^ about 92% of the time at the Born-level 
and about 84% of the time at NLO. These efficiencies are slightly smaller than the b quark 
identification efficiencies shown in Fig. due to the finite widths of the W boson and the 
top quark. 

Fi gure IT31 shows the invariant mass of the reconstructed top quark, comparing the different 
methods to identify the b quark and the W boson from the top quark decay. The top 
quark width is larger than it would be at parton level even though no kinematic smearing 
was applied in this study to mimic the detector effect. This is the result of using the 
reconstructed kinematics of the W boson (in particular the neutrino ^-momentum) rather 
than parton level information. Furthermore, since the same reconstructed W boson was 
used for all curves in Fig. E[ differences among the individual curves are solely due to the 
jet choice used to reconstruct the top quark. Finally, because parton level information is 
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FIG. 15: Invariant mass of the reconstructed W and a jet object. This jet object is: the jet 
containing the b quark from the top decay, using parton information (solid hne), the leading b- 
tagged jet (dashed line), the best jet (dotted line), or the leading jet (dot-dashed line). The W is 
reconstructed either using the standard neutrino p"^ constraint (a) or a top mass constraint (b). 

used for the b jet curve, it functions as a reference and upper limit for the other methods. 

Figure ITM a) shows the top mass reconstructed with the W boson using the standard 
method of choosing the smaller As expected, using the leading jet gives the worst 

performance, and using the leading 6-tagged is a better choice. Given this choice for the 
neutrino, however, the best jet algorithm looks best because for those cases where the 
neutrino has been misreconstructed, the best jet algorithm sometimes chooses the wrong jet 
and thus it falsely appears to give well reconstructed top mass because the algorithm chooses 
an invariant mass close to 178 GeV. Figure ITHT b) shows the top mass reconstructed with the 
W boson reconstructed from a top mass constraint. The overall height of the mass peak is 
higher in Fig. ITH b) than in Fig. ITHT a). indicating that this method is able to reconstruct the 
W boson and b jet correctly more often. It can also be seen that for this choice of W boson 
reconstruction, we can properly reconstruct the true top quark with the leading 6-tagged jet 
with very high efficiency. Using the leading 6-tagged jet together with a top mass constraint 
to reconstruct the W boson gives the best reconstructed top quark and is the most efficient 
event reconstruction algorithm for t-channel single top quark events. We will henceforth 
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FIG. 16: Transverse momentum of the reconstructed top quark, comparing Born-level to 0{as) 
corrections(a) and the individual 0{as) contributions (b). 
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FIG. 17: Rapidity of the reconstructed top quark, comparing Born-level to 0{as) corrections (a) 
and the individual 0{as) contributions (b). 

focus on this algorithm. 

FigurelTBlshows the transverse momentum of the top quark reconstructed from the leading 
6-tagged jet together with the top mass constraint to reconstruct the W boson, hence, to 
determine p'^. The LIGHT and HEAVY corrections tend to shift the top quark to higher 
values, while the TDEC correction lowers it as expected. The different 0{as) corrections 
have only a small effect on the top quark rapidity shown in Fig. El similar to the b jet 
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above. The HEAVY correction tends to shift the top quark to negative rapidities, while the 
LIGHT correction tends to shift it to positive rapidities, resulting in a net effect of only a 
small change in rapidity. 



C. Top Quark Polarization 

Having identified the 6 jet from the top quark decay and the W boson, we can now study 
correlations expected from event kinematics. As noted in the previous section, we shall 
focus on the study of top quark polarization using the leading 6-tagged jet and the top mass 
constraint for the W boson, because it reconstructs the final state correctly with the highest 
efficiency. 

In the SM, the top quark produced in single top quark events is highly polarized, and 
this polarization can in principle be measured. The top quark is by far the heaviest known 
fermion and the only known fermion with a mass at the electroweak symmetry-breaking 
scale. Thus, it is hoped that a detailed study of top quark coupling to other particles will 
be of great utility in determining if the SM mechanism for electroweak symmetry-breaking 
is the correct one, or if new physics is responsible Angular correlations among the 

decay products of polarized top quarks provide a useful handle on these couplings. It has 
been pointed out that among the decay products of the top quark, the charged lepton 
is maximally correlated with the top quark spin We can thus obtain the most 

distinctive distribution by plotting the angle between the charged lepton and the spin axis 
of the top quark in its rest frame. 

Although the top quark is produced via the left-handed charged current, there is no 
reason to believe that the helicity basis will give the best description of the top quark spin. 
Choosing an appropriate basis could maximize spin correlation effects. Three polarization 
bases have been studied in the literature for the t-channel process, and they differ by the 



reference frame used to define t 



le polarization: the helicity basis, the beamline basis, and 



the so-called "spectator" basis All three work in the top quark rest frame, but they 
have different reference axis for the top quark spin, cf. Fig. ^| In the more commonly 
used helicity basis, the top quark spin is measured along the top quark direction of motion 
in the center of mass frame which is chosen as the frame of the (reconstructed top quark, 
spectator jet) system after event reconstruction. In the beamline basis, the top quark spin 
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is measured along the incoming proton direction. In the spectator basis we can maximize 
spin correlations by taking advantage of the fact that the top quark produced through the 
t-channel single top processes is almost 100% polarized along the direction of the spectator 
quark. In the discussion below, we will examine the polarization of single top quark events 
in these three bases. 

In helicity basis, the cm. frame needs to be reconstructed in order to define the top 
quark momentum. Due to additional jet radiation, the determination of the cm. frame at 
NLO is more complicated than at the Born-level. The additional radiation will also blur the 
spin correlation, therefore, choosing the appropriate frame will reduce this effect. In this 
study, two options for reconstructing the cm. frame are investigated: 

1. tg(j)-frame: the cm. frame of the incoming partons. This is the rest frame of all the 
final state objects (reconstructed top quark and all others jets). In exclusive two-jet 
events, this frame is the same as the cm. frame at the Born-level, i.e. reconstructed 
from summing over momentum of the top quark and spectator jet. In exclusive three- 
jet events, this frame is reconstructed by summing over the 4-momenta of top quark, 
spectator jet, and the third-jet from our parton level calculation. 

2. tg-frame: the cm. frame of the top quark and spectator jet. In this case, even in 
exclusive three-jet events, the reference frame is constructed by summing over only 
the 4-momenta of the top quark and spectator jet. Note that this differs from the 
tg(j)-frame only in exclusive three-jet events. 

As shown in Table IIIII and discussed below, the degree of top polarization is larger in the 
tg-frame than in the tg(j)-frame. Therefore, in the figures below we only display the top 
quark polarization in the tg-frame. 

In the helicity basis, the polarization of the top quark is examined as the angular distri- 
bution (cosOhei) of the lepton in the cm. frame relative to the moving direction of the top 
quark in the same frame. The angular correlation in this frame is given by 

\Pt\\Pe*\ 

where p/ is the charged lepton three-momentum defined in the rest frame of the top quark, 
whose three momentum is denoted as pt, which is in turn defined in the cm. frame. For a 
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FIG. 18: Illustration of the three choices for the top quark spin basis. The circle denotes the top 
quark rest frame and the blue arrows denote the top quark spin direction. 




FIG. 19: Top quark polarization in the helicity basis using the full parton information (a) and 
after event reconstruction with selection cuts (b), comparing Born-level to 0{as) corrections. The 
Born-level and NLO curves have been normalized to the same area. 

left-handed top quark, the angular correlation of the lepton is given by (1 — cos 6'/jei)/2, 
and for a right-handed top quark, it is (1 + cos6hei)/2. Figure ITW a) shows that this linear 
relationship for cos Oh^i is indeed a valid description for t-channel single top quark events at 
the parton level. The figure also shows that the top quark is not completely polarized in 
the helicity basis, and that this polarization is only slightly weakened when including 0{as) 
corrections. Figure El(b) shows that this weakening is amplified after event reconstruction, 
where the effect of the lepton-jet separation cut can also be seen, as a drop-off of the cos 6hei 
distribution close to a value of —1. 

In the "spectator" basis, the relevant angular correlation for the t-channel process is 
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FIG. 20: Top quark polarization in the spectator basis using the full parton information (a) and 
after event reconstruction with selection cuts (b), comparing Born-level to 0(qs) corrections. The 
Born-level and NLO curves have been normalized to the same area. 



cos 9 spec, defined as 



cos 6' 
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\Psl 



spec 



(5) 



where Pspec is the spectator jet three-momentum in the top quark rest frame and is the 
lepton three-momentum in the top quark rest frame. Although this basis picks the wrong 
spin axis direction for the db and bd initial states, it is correct most of the time at the Tevatron 
collider. This is because the Tevatron is a pp collider which means that the production rate 
of pp{ub, bu) dt is much larger than pp{db) — > ut. The large slope found in Fig. QUI shows 
that the spectator basis indeed results in a large degree of correlation for t-channel single 
top quark events at the parton level. The slope shown in Fig. I^I¥ a) is opposite to that in 
Fig. irW a) due to the fact that the spin quantization axis points in opposite directions for 
the two basis, cf. Fig. ^1 The degree of polarization of the top quark is larger than that in 
the helicity basis, and the 0{as) corrections blur the spin correlation effects only slightly, 
both at parton level and after event reconstruction. The reconstructed cos Ogpec distribution 
again shows a drop-off due to the lepton-jet separation cut, in this case at high cos 

In the "beamline" basis, the relevant angular correlation for the t-channel process is 
cos ^beam, defined as 

Pp* ■ Pi* 



cos 9, 



beam 



\Pp*\\Pi* 



(6) 
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FIG. 21: Top quark polarization in the beamline basis using the full parton information (a) and 
after event reconstruction with selection cuts (b), comparing Born-level to 0(qs) corrections. The 
Born-level and NLO curves have been normalized to the same area. 

where p* is the proton three-momentum in the top quark rest frame and p/ is the lepton 
three-momentum in the top quark rest frame. In this analysis, we orient the coordinate 
system such that protons travel in the positive z direction and anti-protons travel in the 
negative z direction. For a top quark polarized along the proton moving direction, the 
angular distribution of the lepton is (1 + cos 6 beam) f^, while for a top quark polarized 
along the anti-proton moving direction it is (1 — cos 6'6eam)/2. Figure shows that this linear 
relationship for cos Oheam is a valid description for t-channel single top quark events at the 
parton level. However, the top quark is less polarized in the beamline basis at parton level. 
In this case, the 0{as) corrections actually improve the spin correlation effects at parton 
level. After event reconstruction the situation is similar as before, the spin correlation is 
further reduced and shows a drop-off close to 1. 

To better quantify the change in polarization, it is useful to define the degree of polar- 
ization V of the top quark. This is given as the ratio 

^ N_-N+ 
V = 

iV_ + iV+ ' 

where (A^+) is the number of left-handed (right-handed) polarized top quarks in the 
helicity basis. Similarly, in the spectator (beamline) basis, (A^+) is the number of top 
quarks with polarization against (along) the direction of the spectator jet (proton) three 
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momentum in the top quark rest frame. The angular distribution is then given by j32| 

1 da _ N_ 1 + cos 6^ N+ 1 - cos 6 
ad{cose) ^ N^ + N+ 2 ^ iV_ + iV+ 2 

= i(l + Dcos^,). 

Simple algebra leads to the following identity: 

V = -3 [\-^dx, (7) 
j-i adx 

where — — is the normalized differential cross section as a function of the polar angle x. 
adx 

Here, x denotes cos 9hei in the helicity basis, etc. Based on the degree of polarization V, we 
can easily get the spin fractions J^± as: 

^ l + V 



N^ + N+ 2 ' 



N^ + N+ 2 

Note that JF_(JF_^) is the fraction of left-handed (right-handed) polarized top quarks in the 
helicity basis, etc. 

We can also define the asymmetry A of the distribution as 

^ ^ /°i dajcos 6) - dajcos 6) 

J\ da {cos 9) + /o da {cos 9)' ^ ' 

It is easy to check that without imposing any kinematic cuts, D = 2A. Furthermore, the 
ratio of top quarks with spin along the basis direction will be = 0.5 — A when no cuts are 
applied. However, when cuts are imposed, these two relations break down. Table UTTl shows 
that the relationship V = 2A indeed holds at parton level (within rounding errors) and is 
still approximately true at 0{as)- 

In Table UTTl we present our results for inclusive two-jet events at the parton level before 
selection cuts and after the loose set of cuts and event reconstruction (cf. Sec. IIVB|) . The 
result for exclusive three-jet events is shown in Table ITVl For completeness, we also include 
the same study for the s-channel single top process in the Appendix. 

• We note that in the helicity basis, the degree of top quark polarization is larger in the 
tg-frame than in the tg(j)-frame (the usual cm. frame of the incoming partons) at 
the parton level. The degree of top quark polarization is 94% in the tg-frame and only 
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V 


T 


A 


LO 


NLO 


LO 


NLO 


LO 


NLO 


Helicitv basis: 


Parton level (tofil-frame) 


0.96 


0.74 


0.98 


0.87 


0.48 


0.37 




Parton level (ig-frame) 


0.96 


0.94 


0.98 


0.97 


0.48 


0.47 




Reconstructed events (tg(j)-frame) 


0.84 


0.73 


0.92 


0.86 


0.46 


0.41 




Reconstructed events (tf^-frame) 


0.84 


0.75 


0.92 


0.88 


0.46 


0.42 


Spectator basis: 


Parton level 


-0.96 


-0.94 


0.98 


0.98 


-0.48 


-0.47 




Reconstructed events 


-0.85 


-0.77 


0.93 


0.89 


-0.46 


-0.42 


Beamline basis: 


Parton level 


-0.34 


-0.38 


0.67 


0.69 


-0.17 


-0.19 




Reconstructed events 


-0.30 


-0.32 


0.65 


0.66 


-0.17 


-0.20 



TABLE IIL Degree of polarization "D, polarization fraction JF, and asymmetry A for inclusive 
two-jet single top quark events, at the parton level before cuts and after selection cuts and event 
reconstruction, in the t-channel process. Here, J- corresponds to T- in the helicity basis for left- 
handed top quarks and to Tj^ in the spectator and beamline bases for top quarks with polarization 
along the direction of the spectator-jet and proton three momentum, respectively. Also, the tq{j)- 
frame in the helicity basis denotes the cm. frame of the incoming partons, while the tq-hame 
denotes the rest frame of the top quark and spectator jet. 

74% in tg(j)-frame. This is due to the fact the degree of top quark polarization in 
inclusive two-jet events is a mixture of contributions from both exclusive two-jet events 
and exclusive three-jet events. Table HVl shows that for exclusive three-jet events, the 
degree of top quark polarization is larger in the tg-frame than in the tg(j)-frame. This 
reduction in the degree of polarization for the tg(j)-frame is due to events in which 
the additional jet is produced before the exchange of the virtual W boson. After event 
reconstruction, the two frames give almost the same degree of top quark polarization. 

• We find that at the parton level stage, there is very little difference between the 
helicity basis (using the tg-frame) and the spectator basis, and that both of them 
give significantly better polarization than the beamline basis both at Born-level and 
NLO. The top quark is almost completely polarized in the helicity and spectator bases, 
and the 0{as) corrections only degrade that picture slightly. The similarity between 
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these two bases is due to the fact that the degree of polarization of the top quark is 
dominated by exclusive two-jet events for which these two bases are equivalent, cf. 
the Born-level results in Table 11111 After event reconstruction for the inclusive two- 
jet sample, the degree of polarization is reduced as expected in both the helicity and 
spectator bases. 

• In the beamline basis, the polarization actually increases when going from Born-level to 
NLO, but it is still much lower than in the other two bases. After event reconstruction, 
the degree of polarization is also reduced. 

• In the exclusive three-jet sample, the degree of polarization is further reduced because 
the third-jet affects the kinematics of either spectator jet or the top quark. The helicity 
basis with tg-frame gives almost the same degree of polarization as the spectator basis. 

Our study shows that the helicity basis (using the tg-frame) and the spectator basis are 
equally good to study the top quark polarization. Unlike the s-channel process in which the 
VT-boson is not perfectly reconstructed in the best-jet algorithm and thus the polarization 
measurement was significantly degraded after event reconstruction, using the leading b- 
tagged jet and the top mass constraint gives excellent final state reconstruction in the t- 
channel process, and the degree of top quark polarization is only somewhat degraded after 
event reconstruction. 

D. Distributions for Three-jet Events 

As shown in Fig. IHl a large fraction of the events passing the loose selection cuts contain 
three jets. In this section we focus on the properties of these three-jet events and the 
additional jet. 

In events containing two untagged jets, we can use the pt of the untagged jet to pick up 
the spectator jet. When the b jet from the initial state gluon splitting in the VT-gluon fusion 
process is mistagged, it will also contribute an untagged jet to the event. For simplicity, we 
will assume fully efficient, perfect 6-tagging here and consider events with single 6-tags and 
double 6-tags separately. 

From Fig. IHl it is clear that the jet multiplicity at NLO depends strongly on the jet Pt 
cut. Figure |221 shows that it also depends strongly on the jet pseudo-rapidity cut. The 
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V 




A 


Helicity basis: 


Parton level (tq'(j)-frame) 


0.65 


0.83 


0.35 




Parton level (tg-frame) 


0.78 


0.89 


0.43 




Reconstructed events (tg(j)-frame) 


0.63 


0.81 


0.34 




Reconstructed events (tg-frame) 


0.70 


0.85 


0.38 


Spectator basis: 


Parton level 


-0.78 


0.89 


-0.42 




Reconstructed events 


-0.70 


0.85 


-0.38 


Beamline basis: 


Parton level 


-0.26 


0.63 


-0.16 




Reconstructed events 


-0.27 


0.63 


-0.16 



TABLE IV: Degree of polarization T>, polarization fraction J^, and asymmetry A for exclusive 
three-jet single top quark events, at parton level and after event reconstruction, in the t-channel 
process. Here, corresponds to J-'-, in the helicity basis for left-handed top quarks and to J^^ in the 
spectator and beamline bases for top quarks with polarization along the direction of the spectator- 
jet and proton three momentum, respectively. The tq{j) frame in the helicity basis denotes the 
cm. frame of the incoming partons, while the tq frame denotes the rest frame of the reconstructed 
top quark and light quark. 




FIG. 22: Inclusive cross section and fraction of three-jet events at NLO for varying jet pseudo- 
rapidity cuts, after the loose selection cuts. Shown is the total cross section as a function of the 
jet pseudo-rapidity cut (a) and the fraction of two-jet and three-jet events as a function of the jet 
pseudo-rapidity cut (b). 
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(a) (b) 

FIG. 23: Fraction of events with one or two 6-tagged jets for the inclusive sample and the exclusive 
three-jet sample as a function of the pseudo-rapidity cut on the 6-tagged jets (a), and fraction of 
events with two tagged jets as a function of the pT threshold on the 6-tagged jets for three different 
pseudo-rapidity cuts (b), after the loose selection cuts. 

dependence of the total cross section on the jet pseudo-rapidity cut is different between the 
Born-level and NLO, mostly as a result of the third jet. Figure IHl also shows that not only 
the cross section but also the jet multiplicity depends strongly on the jet pseudo-rapidity 
cut. Only for jet pseudo-rapidity cuts above 3 are cross section and jet multiplicities stable. 

The fraction of 6-tagged jets also depends strongly on the jet pt cut, as shown in Fig. [T21 
There is a large fraction of events in which the third jet comes from the b quark, especially 
for low 6-tag px thresholds. This can also be seen in Fig. ESI which shows the dependence 
of the fraction of events with 1 b-tag and 2 6-tags as a function of the cut on the 6-tagged 
jet Pt and rj. As before, the fraction of events with 2 6-tagged jets is stable only for pseudo- 
rapidity thresholds above 3. At that point, about one quarter of the inclusive events contain 
two 6-tagged jets, and half of the exclusive three-jet events contain two 6-tagged jets. In 
the following analysis, we shall require at least one 6-tagged jet in the event and do not 
distinguish the identity of the third jet, unless specified otherwise. 
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1. Kinematic Distribution of the Extra Jet 

Initial- and final-state emission of additional gluons occurs before the top quark goes on 
shell and can thus be considered as "production-stage emission" , while decay-stage emission 
occurs only after the top quark goes on shell. In principle, an event with an extra jet 
can thus be classified as production-stage or decay-stage by looking at the invariant mass 
of the decay products. In pro duct ion- stage emission events, the W boson and h quark 
momenta will combine to give the top quark momentum. In decay-stage emission events, 
the gluon momentum must also be included to reconstruct the top quark momentum. This 
interpretation is exact at the NLO parton level in the narrow width approximation. Finite 
top width effects can blur the above classification due to interference between production- 
and decay-stage emission. This classification is nevertheless still useful in our case because 
the top width of 1.5 GeV is small compared to the hard gluon Ef cut imposed in the MC 
calculations. It should be kept in mind that in an actual experiment, the production-decay 
distinction is further blurred by the experimental jet energy resolution and ambiguities 
associated with properly assigning partons to jets. 

Figure |^ shows the transverse momentum distribution as well as the pseudo-rapidity 
distribution of the third jet in three-jet events. This jet corresponds to the gluon in about 
70% of the events after the loose set of cuts. The transverse momentum distribution of 
the third jet for those events where that jet corresponds to the h quark (from VT-gluon 
fusion subprocess) can be seen in Fig. El It comprises about 80% of the HEAVY correction 
after imposing the loose selection cuts, which dominates over LIGHT and TDEC radiative 
corrections. 

As expected, the Et distribution is steeply falling for all contributions, but it extends 
to much higher values for HEAVY emission. The smaller values of Et to which TDEC 
emission is constrained are again a consequence of the available phase space from the top 
quark decay. Recall that the top quarks are produced with relatively modest transverse 
momentum (cf. Fig. fTH|) . so that gluons from the decay do not receive much of a boost. 
Note also that an increase in the pt cut on the jet would result in a further reduction in 
relative size of the decay contribution compared to production. Figure |21] also shows the 
distribution in pseudo-rapidity of the extra jet. The third jet radiated from the LIGHT 
and HEAVY quark lines has a relative larger magnitude in its (non-symmetric) rapidity. 
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FIG. 24: Transverse momentum (a) and pseudo-rapidity (b) of the third jet after selection cuts for 
the various 0(a,) contributions. 
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FIG. 25: Separation between the third jet and the tagged jet (a) and between the third jet and 
the untagged jet (b) after selection cuts for the various 0{as) corrections. 

as compared to the more central TDEC emission. This is consistent with our intuition 
that decay-stage radiation, dominated by the gluon radiated from the bottom quark which 
tends to appear in the central pseudo-rapidity region, is also likely to be produced centrally. 
However, this TDEC contribution is small and the HEAVY radiation dominates even in the 
central region. 

This tendency of decay-stage radiation to be associated with the final-state b quark might 
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lead one to expect that if the extra jet is "near" the b jet it should be included in the top 
quark reconstruction, and if it is not then it should be excluded. Figure 123 which shows 
the angular separation AR between the extra jet and the leading 6-tagged-jet as well as 
the leading untagged jet, confirms that the decay-stage radiation peaks close to the leading 
6-tagged jet, and pro duct ion- stage radiation peaks farther away. Figure O^T a) clearly shows 
that the decay contribution dominates in the low AR region. This is different from the 
s-channel single top process, in which the production contribution dominates over decay 
emission even in the small AR region, cf. Fig. 17 of Ref. jl^. A higher cut on the jet 
would make this situation worse because it would increase the relative size of the production 
emission. 

Nevertheless, the figure suggests that it might be possible to further improve the top 
quark reconstruction for exclusive three-jet events. In this paper, we have been using the 
leading 6-tagged jet to identify the b quark from the top quark decay, and constraints on 
and Mw to obtain the correct p^- With this procedure the final state can be reconstructed 
accurately about 84 % of the time in inclusive two-jet events, cf. Fig. EI Part of the 
efficiency loss is due to the gluon radiated in TDEC emission, and Fig. 123 indicates that 
some of this loss could be reclaimed by including the third jet in the top quark reconstruction 
if it is close to the 6-tagged jet. The actual cut on AR would need to be tuned to maximize 
the top quark reconstruction efficiency. However, tuning a prescription for dealing with the 
extra jet in t-channel single top quark events is complicated because the effects of multiple 
emission, hadronization, and detector resolutions will affect the result. For simplicity, we 
thus chose not to include the third jet in the top quark reconstruction algorithm used in 
this paper. 

In Fig. I^^ b). the equivalent distribution in AR between the extra jet and the untagged 
jet is also shown, where the LIGHT quark line radiation peaks close to the untagged jet as ex- 
pected from gluon radiation in the final state. Obviously, the radiation of the HEAVY quark 
line is dominated by initial state radiation, therefore, its contribution to Ai?(j3, juntagged) is 
small in the region of small AR. The peak at zero in this distribution corresponds to events 
containing two 6-tagged jets in the the W^-gluon fusion subprocess. 
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FIG. 26: Angular correlation cos 9 between the gluon and the b quark before any selection cuts 
using the full parton information (a) and between the third jet and the best jet after selection cuts 
(b). The solid line shows all 0{as) contributions except for the decay part, while the dashed line 
shows only the 0{as) decay contribution. 

2. Angular Correlation Between the Extra Jet and the Best Jet 

As discussed above, the Ai? separation between the third jet and the leading 6-tagged 
jet can be used to distinguish decay-stage emission from pro duct ion- stage emission. It was 
shown in Ref. Q| that the best-jet algorithm is effective in reconstructing the s-channel 
single top quark event, while retaining the top quark polarization information, but can also 
distinguish the decay-stage gluon radiation from the pro duct ion- stage radiation by studying 
which jets are chosen as the best jets. Similarly, we can also study the correlations between 
the extra jet and the best jet in t-channel single top quark events using the best jet algorithm. 

Figure 1211 (a) shows the angular correlation cos 9 between the radiated parton (gluon or 
h) and the h quark at parton level before cuts. Figure OBI fb) shows the same correlation after 
event reconstruction between the third jet and the best jet. Only events for which the best 
jet algorithm chooses a two-jet system are included in the figure. In this case there is a clear 
separation between production-stage and decay-stage emission, and the best jet algorithm 
can be used to separate the two. This approach thus allows a detailed experimental study 
of the radiation pattern in top quark decays in the t-channel single top quark events. 
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V. CONCLUSIONS 



We have presented a next-to-leading order study of t-channel single top quark events at 
the Tevatron, including 0{as) QCD corrections to both the production and decay of the top 
quark. To obtain an accurate prediction of the inclusive rate of t-channel single top quark 
production, a modified narrow width approximation has been adopted to link the production 
of the top quark with its decays (thus preserving top quark spin information) instead of the 
usual narrow width approximation. The impact of kinematical cuts on the acceptances has 
been studied for several different sets of cuts. We found that the acceptances are sensitive 
to the ARcut we imposed on the jet cone size and the lepton isolation. With the choice of 
Ai?cut = 0.5, we found that the difference between the Born-level and NLO acceptances is 
about 10% for a loose cut set (a) and 3 ~ 6% for a tight cut set (c). The above difference 
becomes significantly larger when changing Ai?cut from 0.5 to 1.0. 

We categorize the 0{as) contributions to the t-channel single top process into three gauge 
invariant sets: the light quark line corrections, the heavy quark line corrections and the top 
quark decay corrections. Keeping track of the different categories facilitates the comparison 
between event generators and exact NLO predictions. 

The 0{as) corrections are small in size and contribute about 5.5% of the inclusive cross 
section at NLO. They do however modify the event kinematics and in particular result in a 
large fraction of events containing three reconstructed jets in the final state for the loose set 
of kinematic cuts. The acceptance for t-channel single top quark events and this fraction of 
3-jet events depend strongly on the pseudo-rapidity cut on the jets. Although the radiative 
corrections to the inclusive t-channel single top production rate are small, they affect the 
shape of some of the important kinematic distributions that can be used experimentally to 
separate the t-channel single top signal from the various backgrounds, such as the pseudo- 
rapidity distribution of the spectator jet. We find that the 0{as) LIGHT and HEAVY 
corrections have almost opposite contributions to various pseudo-rapidity distributions, due 
to the difference in the parton distribution functions between the valence quarks and sea 
quarks. The former shifts the spectator jet to even higher pseudo-rapidities, while the 
later shifts it to more central pseudo-rapidity regions. The summed contributions cause the 
spectator jet to be even more forward which will change the prediction of the acceptance 
for t-channel single top quark events. Also, a large fraction of three jet events contain two 
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fe-jets due to the collinear enhancement in the W + g fusion process. This imphes that higher 
order corrections need to be calculated in order to correctly predict the behavior of the 6-jet 
in the small region. 

In order to study top quark properties such as the top quark polarization, induced from 
the effective t-b-W couplings, we need to reconstruct the top quark by combining the re- 
constructed W boson with the 6-tagged jet. Most of the t-channel single top quark events 
contain only one 6-jet in the final state. Thus, we can use the leading 6-tagged jet algorithm 
to identify the 6-jet in the final state, rather than the best-jet algorithm which is more 
appropriate for the s-channel. We found that the leading 6-tagged jet algorithm effectively 
picks up the correct 6-jet in the event. Since this algorithm makes no assumptions about 
the other particles in the event, we can use it to also determine the longitudinal momentum 
of the neutrino (p^) accurately through a top quark mass constraint. Information about 
the reconstructed final state can be used to explore correlations between different objects 
in the event. After reconstruction of the top quark event, we study spin correlations in the 
top quark decay in three different bases: the helicity basis, the spectator basis, and the 
beamline basis. We find that the degree of top polarization is very large, especially in the 
helicity and spectator bases. This is true even after event reconstruction because we are 
able to reconstruct the top quark final state almost perfectly. As one expects, the degree of 
top polarization is reduced slightly by the 0{as) corrections, both at parton level and after 
event reconstruction with the loose selection cuts. We also note that using the tg-frame 
(the rest frame of the reconstructed top quark and spectator jet) in the helicity basis gives 
almost the same degree of polarization as in the spectator basis. 

Finally, we point out that the above conclusion does not strongly depend on the choice of 



jet algorithm. We have checked that at the parton level using the Durham algorithm 
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3J] leads to similar conclusion on the relative importance of the production- and decay-stage 



gluon emission. 
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APPENDIX A: TOP QUARK POLARIZATION IN THE S-CHANNEL PROCESS 



Above we presented the top quark spin correlations in the hehcity basis in the t-channel 
using two different definitions for the cm. frame. In our previous paper on the s-channel, 
the degree of top quark polarization in the helicity basis was calculated only in the tb frame 
(see Sec. IV of Ref. and given only for inclusive 2-jet events. That frame corresponds 
to the cm. frame of the incoming partons strictly speaking only for events without initial 
or final state gluon radiation. For completeness, we also present here the results for the top 
quark polarization in the s-channel single top process for two different cm. frame definitions. 
Similar to the t-channel single top process, we define two reference frames to study the degree 
of top quark polarization at the parton level in the helicity basis for s-channel single top 
quark events. 

1. t6(j)-frame: the cm. frame of the incoming partons. This is the rest frame of all 
the final state objects (reconstructed top quark and all others jets). In exclusive two- 
jet events, this frame is the same as that at the Born-level, i.e. reconstructed from 
summing over momentum of the top quark and non-best-jet. In exclusive three-jet 
events, this frame is reconstructed by summing over the 4-momenta of top quark, 
non-best-jet, and the third-jet from the parton level calculation. 

2. t6-frame: the cm. frame of the top quark and non-best-jet. In this case, even in the 
exclusive three-jet events, the reference frame is constructed by summing over only 
the 4-momenta of the top quark and non-best-jet. Note that this differs from the 
i:6(j)-frame only in exclusive three-jet events. 

We present the resulting top quark polarization in the helicity and optimal bases in 
Table |V] for inclusive 2-jet events and in Table IVII for exclusive 3-jet events. Similar to 
the t-channel, the degree of polarization is larger in the helicity basis in the tb frame. The 
polarization is reduced when including the possible third jet to reconstruct the cm. frame in 
particular for those events where that jet is produced in the initial state, before the virtual 
W boson is created. Table IVTl also shows that the degree of polarization is generally smaller 
for 3-jet events, which is again similar to the t-channel. 
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V 




A 






LO 


NLO 


LO 


NLO 


LO 


NLO 


Helicity basis: 


Parton level (t6(j)-frame) 


0.63 


0.54 


0.82 


0.77 


0.32 


0.27 




Parton level (t6-frame) 


0.63 


0.58 


0.82 


0.79 


0.32 


0.29 




Reconstructed events (t6(j)-frame) 


0.46 


0.37 


0.73 


0.68 


0.21 


0.21 




Reconstructed events (t6-frame) 


0.46 


0.37 


0.73 


0.68 


0.26 


0.21 


Beamline basis: 


Parton level 


-0.96 


-0.92 


0.98 


0.96 


-0.48 


-0.46 




Reconstructed events 


-0.48 


-0.42 


0.74 


0.71 


-0.24 


-0.21 



TABLE V: Degree of polarization polarization fraction J^, and asymmetry A for inclusive 
two-jet single top quark events, at the parton level and after event reconstruction, in the s-channel 
process. Here, J- corresponds to T- in the helicity basis for left-handed top quarks and to in the 
beamline basis for top quarks with polarization along the direction of anti-proton three momentum, 
respectively. The thg frame in the helicity basis denotes the cm. frame of the incoming partons 
while th frame denotes the rest frame of the reconstructed top quark and b quark. 





V 




A 


Helicity basis: 


Parton level {thg frame) 


0.45 


0.72 


0.25 




Parton level {tb frame) 


0.49 


0.74 


0.27 




Reconstructed events {thg frame) 


0.36 


0.68 


0.21 




Reconstructed events {th frame) 


0.37 


0.68 


0.21 


Beamline basis: 


Parton level 


-0.81 


0.91 


-0.45 




Reconstructed events 


-0.38 


0.69 


-0.19 



TABLE VI: Degree of polarization P, polarization fraction and asymmetry A for exclusive 
three-jet single top quark events, at parton level and after event reconstruction, in the s-channel 
process. Here, J- corresponds to J-^ in the helicity basis for left-handed top quarks and to J-j^ in the 
beamline bases for top quarks with polarization along the direction of anti-proton three momentum, 
respectively. The thg frame in the helicity basis denotes the cm. frame of the incoming partons 
while the tb frame denotes the rest frame of the reconstructed top quark and b quark. 
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